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Abstract. We use the very deep and homogeneous I-band selected dataset of the FORS Deep Field (FDF) to trace 
^ ' the evolution of the luminosity function over the redshift range 0.5 < z < 5.0. We show that the FDF 1-band selec- 

l/^ I tion down to Iab = 26.8 misses of the order of 10 % of the galaxies that would be detected in a K-band selected 

CO . survey with magnitude hmit Kab ~ 26.3 (hke FIRES). Photometric redshifts for 5558 galaxies are estimated 

' based on the photometry in 9 filters (U, B, Gunn g, R, 1, SDSS z, J, K and a special filter centered at 834 nm). 

, A comparison with 362 spectroscopic redshifts shows that the achieved accuracy of the photometric redshifts is 

Az/{z3pec + 1) < 0.03 with only ~ 1% outliers. This allows us to derive luminosity functions with a refiabihty 
similar to spectroscopic surveys. In addition, the luminosity functions can be traced to objects of lower luminosity 
which generally are not accessible to spectroscopy. We investigate the evolution of the luminosity functions eval- 
• uated in the restframe UV (1500 A and 2800 A), u', B, and g' bands. Comparison with results from the hterature 

shows the reliabihty of the derived luminosity functions. Out to redshifts of 2: ~ 2.5 the data are consistent with 
Q . a slope of the luminosity function approximately constant with redshift, at a value of —1.07 ± 0.04 in the UV 

$H ' (1500 A , 2800 A) as well as u', and —1.25 ±0.03 in the blue (g', B). We do not see evidence for a very steep slope 

I (a < —1.6) in the UV at (z) ~ 3.0 and (z) ~ 4.0 favoured by other authors. There may be a tendency for the faint- 

ed ■ end slope to become shallower with increasing redshift but the effect is marginal. We find a brightening of M* and 

a decrease of (p* with redshift for all analyzed wavelengths. The effect is systematic and much stronger than what 
can be expected to be caused by cosmic variance seen in the FDF. The evolution of M* and (p* from z = to 2 = 5 
is well described by the simple approximations M*{z) = Mq + aln(l + z) and 4>*{z) — (f>o{l + z'f for M* and 0*. 
5— ( ' The evolution is very pronounced at shorter wavelengths (a = —2.19, and h — —1.76 for 1500 A rest wavelength) 

and decreases systematically with increasing wavelength, but is also clearly visible at the longest wavelength 
investigated here (a = —1.08, and b = —1.29 for g'). Finally we show a comparison with semi-analytical galaxy 
formation models. 

Key words. Galaxies; luminosity function - Galaxies: fundamental parameters - Galaxies: high-redshift - Galaxies: 
distances and redshifts - Galaxies: evolution 
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1. Introduction 

Observational constraints on galaxy formation have im- 
proved significantly over the last years and it has become 
possible to study the evolution of global galaxy proper- 
ties up to very high redshifts. A crucial step to probe 
the properti es of galaxies up to t he hig hest redshifts was 
the w ork of ISteidel fc HamiltonI l)l993l) and ISteidel et alJ 
l|l996^ who used color selection to discriminate between 
low redshift and high redshift galaxies. Although the 
Lyman-break technique is very efficient in selecting high 



2 



Gabasch et al.: The evolution of the luminosity functions in the FDF: I. The blue bands 



redshift galaxies fsee lBlaizot et alJboOSl for a detailed dis- 
cussion) with a minimum of photometric data, it has the 
disadvantage that it does not sample galaxies homoge- 
neously in redshift space and may select against certain 
types of objects. With the advent of deep multi-band 
photometric su r veys (Hubble Deep Field North (HDFN 



ICole et all hOOd iMenci et al.1. I2002L 12003^ and simula- 



Williams et al. 
Arnouts et al.l 



Il99fil) . NTT SUSI deep Field (NDF 
1999t) Hubble Deep Fie l d Sou th (HDFS 



WilUams et alJ 12000; "C asertano et all l200rf >. Chandra 
Deep Field South fCDFS: lArnouts et a,lJl2nnHl. W illia.m 
Herschel Deep Fiel d (WHDF: iMcCracken et all l2000t 
Metcal fe et all l200ll) Subaru Deep Fie ld/Survey (SDF; 
Maihara et alJl200lHOuchi et al.ll2003a|) . 

The G reat Observator ie s Or igins Deep Survey 
(GOODS; iGiavalisco et all hOO^ ) the photometric 
redshift technique (essentially a generalization of the 
drop-out technique) has increasingly been used to 
identify high-redshift galaxies. Several methods have 
been descr ibed in the li te ratur e to de ri ve photometric 
redshift,s JrJ ITflBi Ik^ ITflsi iRnirmer et al. , 
1995 iFornandez-Soto et a,l.L M \BeniteA l200r : 
Le Borgne fc Rocca-Volmerangel l2002t iFirth et al. . 



20031) . 



Based on either spectroscopic redshifts, drop-out 
techniques, or photometric redshifts, it has been possible 
to derive luminosity functions at different redshifts in 



19991 ICowie et al. 
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; .Adelberger & SteideL .200a 


Cohen et all l200( ; 


ISullivan et all l2000l lOuchi et al. 


2001 


|Poli et al.' '2001; 


"Wilson et al.' '2002; 'Wolf et al 
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iRowan-Robinson^ 




2003; .Kashikawa et al 2003 



Onchi el 



bands llLillv et all 11995; 'Heyl ct al.', 'l9 97t iLin et al 
1997t lSawicki et all fl997.;. Small et al 1997HZucca et al 
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all il999t lMa,rmoni et a,].L 



ii^oyeaav eL aiA , 
Fried et all '2001: Cross & Driver, 2002, .... 
Marinoni et al.. 2002: .Norberg et al , 2002: BcU ct al 
20031 Ide Lapparent et all 120031: iLiske et all I2OO' 
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'im ct al.", 2002 ; 



PoU et al.U2003HPerez-Gonzalez et al.ll2003 'l. Within the 
uncertainties given by IMF and dust content, the flux 
in the UV al l ows t o trace the star formation rate (SFR; 
iMadau et al.l Il99 8^ in the galaxies, while the optical 
luminosities provide constraint s on more evolved stellar 
populations l)Franx et al.U200^ . 

Lo cally, the 2dF G alaxy Redshift Survey (2dF- 



GRS; IColless et a,].1l200lh . the Sloan Digital Sky Survey 
(SPSS: IStonghton et al I [ 

vey l|,Tarrett et al.L I2OO0I) have provided superb refer- 



hom and the 2MASS 



sur- 



ence points for galaxy l uminosity funct i ons ov er a large 
wavele ngth range ( s ee iNorberg et alJ l2002l for 2dF- 
GRS, iB lanton et al.' '200 1', '2003' for the SDSS, and 
iKochanek et al. 2001: QqIc et al. 2001 for 2MASS). 

In parallel to the observational effort, theoreti- 
cal models have been developed within the frame- 
work of the cold dark matter cos mology. Most no- 
tably , semi-analyt i c mod e ls (SAMs) feauffmann et all 
1993; 'Cole et all Il 994: 'SomcrviUc & Primacl^, 'l999l 
Kauffmann et allll999;.Poh et al.. .1999.: -Wu et ak .20001 



tions based 


on 


smoothed-particle hydrodynamics 




(IDave et al. 


"1999: 'Weinberg et alJ 


l2002t 


Nagamind 


(SPHj 


iNagamine et al.. 


.2003il have made testable predictions 



Starting with the mass function of dark matter halos and 
their merging history, SAMs use simplified recipes to de- 
scribe the baryonic physics (gas cooling, ph otoionization, 
star f ormation, feedback processes, etc., see lBenson et alJ 
I2OO3I) to derive stellar mass and luminosity functions. 

Ideally, a comparison between observations and models 
should be done with deep multiwavelengths datasets that 
also cover a large area. The dataset has to be sufficiently 
deep in order to be able to derive the faint-end slope of 
the luminosity function. On the other hand, one also needs 
as large an area as possible to overcome cosmic variance 
and to quantify the density of rare bright galaxies, which 
define the cut-off of the luminosity function. 

The FORS Deep Field llHeidt et alll2003l) has a depth 
close to the HDFs but an area of 8 - 10 times the area 
of the HDFN. This depth allows us to detect galaxies at 
z > 2 which would be missed by Lyman-b reak studies 
which us ually reach only Ra r < 25 .5 (see also lFranx et alJ 
I2OO3I and Ivan Dokkum et alJl2003|) . 

Very reliable photometric redshifts are crucial for 
the analysis of the evolution of the luminosity func- 
tions in the FDF. Photometric redshifts have been de- 
ter mined with a templ ate matching algorithm described 
m iBender et af] l|200j) that applies Bayesian statistics 
and uses semi-empirical template spectra matched to 
broad band photometry. We achieved an accuracy of 
Az/{zspec + 1) < 0.03 with only ~ 1% extreme outliers 
(numbers based on a comparison with 362 spectroscopic 
redshifts). Redshifts of galaxies that are several magni- 
tudes fainter than typical spectroscopic limits could be 
determined reliably and thus allowed better constraints 
on the faint-end slope of the luminosity functions. 

In this paper we present the redshift evolution of the 
luminosity function evaluated in the restframe UV-range 
(1500 A, 2800 A), u' (SDSS), B, and g' (SDSS) bands in 
the redshift range 0.5 < z < 5.0. Luminosity functions at 
longer wavelengths as well as the evolution of the luminos- 
ity density and the star formation rate will be presented 
in future papers (Gabasch et al., in preparation). We pro- 
vide a short description of the FDF in Sect. [3 where we 
also present the selection criteria of our galaxies. In Sect. 01 
we investigate possible selection effects due to our purely 
I-band selected catalogue. In Sect.0]we discuss the accu- 
racy of the photometric redshifts as well as the redshift 
distribution of the selected galaxies. In Sect. Eland in the 
Appendix we show luminosity functions at different wave- 
lengths and redshifts. In Sect. El a p arameterization o f the 
redshift evolution of the Schechter (jSchechteil Il976^ pa- 
rameters M* and (j)* is given. We compare our results with 
previous observational results in Sect. [71 and with model 
predictions in Sect. |H1 before we summarize this work in 
Sect. El 
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We use AB magnitudes and adopt a A cosmology In total 78 objects were classified as stars and removed 
throughout the paper with f^M = 0.3, flA — O.l, and from our sample. Our final I-band selected catalogue com- 
Hq = 70kms^^ Mpc^^. prises therefore 5558 objects. 



2. The FORS deep field 

The FORS Deep Field l|Appenzeller et all l2000(l is a 

multi-color photometric and spectroscopic survey of a 
7' X 7' region near the south galactic pole including the 
QSO Q 0103-260 at redshift z = 3.36. The data have been 
taken with FORSl and F0RS2 at the ESO VLT and SofI 
at the NTT. 

The data in the U, B, g, R, I, J and Ks filters were 
reduced and calibrated (inclu ding the correction for galac- 
tic extinction) as described in lHeidt et alJ l)2003j) . The re- 
duction of the images in the z-band and the special filter 
centered at 834 nm follows the same recipe, except for 
additional de-fringing in the z-band. 

The images were stacked with weights to get optimal 
signal to noise for point-like faint objects. The formal 50% 
completeness limits for point sources are 26.5, 27.6, 26.9, 
26.9, 26.8, - 25.5, ~ 25.8, 23.8, 22.6 in U, B, g, R, I, 
834 nm, z, J and Ks, respectively. The seeing varied from 
0.5 arcsec in the I and z band to 1.0 arcsec in the U-band. 
Because the depth of the images decreases towards the 
borders, we limited our analysis to the inner 39.81 arcmin^ 
of our field. The signal-to- noise ratio (S/N) in this 'deep' 
region is at least 90 % of the best S/N in every filter. This 
prevents a possible bias of the photometric redshifts (see 
Sect.EJ due to a not completely homogeneous dataset. 

Object d etection was done in the I-band image using 
SExtractor l|Bertin & Arnout and the catalogue 

for this 'deep' part of the FDF includes 5636 objects. To 
avoid contamination from stars, we rely on three sources 
of information: The star-galaxy classifier of the detection 
software SExtractor, the goodness of fit for galaxy ob- 
jects of the photometric redshift code and, if available, on 
the spectroscopic information. We first exclude all bright 
(J < 22™) starlike objects (SExtractor star galaxy classi- 
fier > 0.95). Then we exclude all objects whose best fit- 
ting stellar spectral energy distribution (SED) - according 
to the photometric redshift code - gives a better match 
to the flux in the different wavebands than any galaxy 
template (2 xltar < X^aiaxy)- These objects are subse- 
quently flagged as star and removed from our catalogue. 
Further inspection of the images confirms, that none of 
these flagged objects are extended. Finally, we reject all 
objects spectroscopically classified as stars. We checked 
the influence of misidentified or missed stars on the lu- 
minosity functions. If stars are fitted by galaxy templates 
their redshifts are mostly very small {z < 0.15, especially 
if they are G and K stars) and, therefore, did not enter 
the analysis. M stars interpreted as galaxies tend to be 
distributed more evenly in redshift space but they do not 
contribute significantly to the number density in any red- 
shift interval. Even if all stars were included as galaxies in 
the sample, they would not affect the derived luminosity 
functions at a noticeable level. 



3. I selection versus K selection 

We use the ultradeep near-infrared I SAAC observat ions 
of the Hubble Deep Field South (,Labbe et all l2003l) for 
a more quantitative analysis of possible selection effects 
between K and I band selected samples. 

In Fig. n we show the Isi4 — Kg versus Kg color- 
magnitude relati on for i^^.^-se le cted o bjects of the HDF- 
b as given by iLabbe et all ((20031) (data were taken 
from : http:/ /www .strw.leidenuniv.nl/~fires/). Following 
iLabbe et al.l l|2003l) , only sources with a minimum of 20% 
of the total exposure time in all bands are included and 
shown as filled symbols. Colors are plotted with la error 
bars. The solid line corresponds to the 50 % completeness 
limiting magnitude of the FDF in the I-band (/ ~ 26.8). 
The figure clearly shows, that although we selected in I, 
we miss only about 10 % of the objects that would have 
been detected in deep K-band images (with a 50 % com- 
pleteness limiting magnitude of Kab ~ 26.3). All objects 
on the left of the solid line would have been detected in the 
I-selected FDF catalogue as well. Therefore we conclude 
that only a small fraction (~ 10 %) of galaxies is missed 
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Fig. 1. /gi4 — Ks versus Ks color-magnitude relatio n for 
JC.s-se lected objec ts of the HD F-S as given bv .Labbc et ^ g,lj 
(|2003l) . Following FLabbe et ai.l l|200^ onlv sources with a 
minimum of 20% of the total exposure time in all bands 
are included and shown as filled symbols. Colors are plot- 
ted with ItT error bars. The solid line corresponds to the 
limiting magnitude of the FDF (1=26.8). Only the objects 
on the right of the solid line are beyond our I-band limit. 
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in deep I-band selected samples relative to deep K-band 
selected samples, provided the I-band images are about 
0.5 AB-magnitudes deeper than the K-band images. Of 
course, this holds only for galaxies at redshift below 6. At 
higher redshifts no signal is detectable in the I-band, due 
to the Lyman break and intervening intergalactic absorp- 
tion. 

Another indication that we are unlikely to miss a large 
population of high redshift red galaxies comes from Fig. 01 
(left panel). Out to redshifts of about 1.5, red galaxies 
define the bright end of the luminosity function. Beyond 
z ~ 1.5 bluer star-forming galaxies take over. Red galaxies 
could still be detected at z > 1.5 but seem to be largely 
absent. In any case, even if we missed a few objects, the 
evolution of luminosity functions that we discuss below 
will not be affected. 

As a side remark we note that also a B-band selected 
FDF catalogue delivers similar conclusions on the evolu- 
tion of the luminosity functions out to redshift ~ 3. Again, 
above this redshift no signal is detectable in the B-band 
due to the Lyman break and intervening intergalactic ab- 
sorption. 

4. Photometric redshifts 

A brief summary of the photometric redshift technique 
used to derive the distances t o the galaxies in the FDF 
can be found in iBender et al.l l)200l|) , a more detailed de- 
scription will be published in a future paper (Bender et al. 
2004). Well determined colors of the objects which implies 
very precise zeropoints in all filters are crucial to derive 
accurate photometric redshifts. Therefore we checked and 
fine-tuned the calibration of our zeropoints by means of 
color-color plots of stars. We compared the colors of FDF 
st ars with the co lors of stellar templates from the library 
of lPicklesI l(l998^ converted to the FORS filter system. In 
general, corrections to the photometric zeropoints of only 
a few hundredth of a magnitude were needed to obtain an 
optimal match to the stars and best results for the pho- 
tometric redshifts. In order to avoid contamination from 
close-by objects, we derived object fluxes for a fixed aper- 
ture of 1.5" (1.5 X seeing) from images which had been 
convolved to the same point spread function. A redshift 
probability function P(z) was then determined for each 
object by matching the object's fluxes to a set of 30 tem- 
plate spectra redshifted between z = and z = 10 and 
covering a wide range of ages and star formation histo- 
ries. As templat e s we u sed ( a) local galaxy templa tes from 
iMannucci etall l|200lj) . and lKinnev et all l)l996l) and (b) 
semi-empirical templates more appropriate for modest to 
high redshift galaxies. The semi-empirical templates were 
constructed by fitting combinations of theore tical spec- 
tral en ergy distributions of diff e rent a ges from iMarasto":] 
ijiggS*) and 'Bruzual & CharloJ l|l993^ with variable red- 
dening ( Kiimey_ct_ah , ,1994) to the observed broad band 
colors of about 100 galaxies in the Hubble Deep Field 
and about 180 galaxies from the FDF with spectroscopic 
redshifts. The remaining 180 galaxies in the FDF with 



spectroscopic redshift were used as an independent con- 
trol sam ple. Lyman for est absorption was parameterized 
following iMadaiJ 1 19951) and references therein. 

In Fig. 121 we compare the photometric and spectro- 
scopic redshifts o f 362 galaxies and QSOs in the FDF (see 
INoU et all [20031 iBohm et alJl2003l for the spectroscopic 
redshifts). The agreement is very good and we have only 
6 outliers with a redshift error larger than Az > 1 among 
362 objects. Three of the outliers are quasars or galaxies 
with a strong power-law AGN component (crosses). The 
others are very blue objects with an almost feature- less 
continuum (triangles). Fig. (left panel) presents the 
distribution for the best fitting template and photometric 
redshifts. Note that to calculate the x^, we have used the 
observational photometric errors and, in addition, have 
assumed that the templates have an intrinsic uncertainty 
of typically 5% in the optical bands and 20% in the in- 
frared bands. The larger errors for the near-IR take into 
account the slightly lower quality of the infrared data if 
compared to the optical. Allowing for this intrinsic uncer- 
tainty turns a discrete set of templates into a template- 
continuum. Observational errors and intrinsic 'errors' were 
added in quadrature. The median value of the reduced 
X^ is below 1.7 and demonstrates that the galaxy tem- 
plates describe the vast majority of galaxies in the FDF 
very well. The right panel of Fig.Oshows the distribution 
of the redshift errors. It is nearly Gaussian and scatters 
around zero with an rms error of Az/{zspec -I- 1) ~ 0.03. 
In Fig. ^ (left panel) , we plot the absolute B-band mag- 
nitudes against the photometric redshifts of the objects. 
Colors from red to blue indicate increasingly bluer spec- 
tral energy distributions. The two lines indicate the 50% 
completeness limit for a red and a blue spectral energy dis- 
tribution corresponding to an I-band limiting magnitude 
of 26.8. The redshift histogram of all objects in the FDF 
is shown in the right panel of Fig. 01 (see also Table 
Most if not all peaks in the distribution are due to real 
clustering in redshift space. From the 362 spectroscopic 
redshifts, we have identified clusters, groups or filaments 
of galaxies with more than 10 identical or almost identi- 
cal redshifts at z 0.22, z = 0.33, z = 0.39, z = 0.45, 
z = 0.77, z — 2.35. Other structures (with only a few 
identical spectroscopic redshifts) are possibly present at 
z = 0.95, z ^ 3.15, and z = 3.4. 



5. Luminosity functions 

5.1. The method 

We compute the absolute magnitudes of our galaxies us- 
ing the I-band selected catalogue as described in Sect. (21 
and the photometric redshifts described in Sect. 01 To de- 
rive the absolute magnitude for a given band we use the 
best fitting SED as determined by the photometric red- 
shift code and convolve it with the appropriate filter func- 
tion. As the SED fits all 9 observed-frame wavebands si- 
multaneously, possible systematic errors which could be 
introduced by using K-corrections applied to a single ob- 
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Fig. 2. Comparison of spectroscopic ((Noll et all l200.'4 iBohm et all 1200,'^ and photometric redshifts for different 
galaxy types and quasars in the FDF (362 objects). 



served magnitude are reduced. Since the photometric red- 
shift code works with 1.5" aperture fluxes, we only need to 
correct to total luminosities by applying an object depen- 
dent scale factor. For this scale factor we used the ratio 
of the I-band aperture flux to the total flux as provided 
by SExtractor (MAG_APER and MAG_AUTO). We have 
chosen the I-band because (a) our I-band data are very 
deep, (b) all objects were detected and selected in the 
I-band, and (c) high redshift galaxies have only poorly 
determined or no flux at shorter wavelengths. This pro- 
cedure may introduce a slight bias, as galaxies are more 
compact or knotty in the rest-frame UV bands (tracing 
HII regions) than at longer wavelengths. However, scal- 
ing factors derived in the deep B-band turned out to be 
similar (for low enough redshifts). 



In a given redshift interval, the luminosity function 
is computed by dividing the number of galaxies in each 
magnitude bin by the volume Vbin of the redshift inter- 
val. To account for the fact that some fainter galaxies 
are not visible in the w' hole survey volume we perform a 
V/Vmax llSchmidtLll968|) correction. Using the best fitting 
SED we calculate the maximum redshift Zmax at which 
the object could have been observed given the magnitude 
limit of our field. We weight each object by Vbin/Vmax 
where Vtm is the volume of our redshift bin enclosed by 
ziow and Zhigh and Vmax is the volume enclosed between 

To derive reliable Schechter parameters we limit our 
analysis of the luminosity function to the bin where the 
V/Vmax begins to contribute at most by a factor of 3 (we 
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Fig. 3. Left panel: Histogram of the reduced for all galaxies in the FDF as obtained for the best fitting template 
and redshift. The dotted vertical line indicates the median reduced x^- R-ight panel: Histogram of the photometric 
redshift errors. The error distribution can be approximated by a Gaussian centered at 0.006 with an rms of 0.028 
(dotted line). 



Table 1. Galaxy distribution in the FDF for the redshift 
intervals used for computing the luminosity function. Note 
that we derive the luminosity function in all redshift bins, 
but exclude the lowest (z < 0.45) and highest redshift 
bin (z > 5.01) from our analysis of the luminosity func- 
tion evolution, since the lowest redshift bin corresponds 
to a too small volume while the z > 5.01 bin suffers from 
incompleteness . 



redshift 




number 


fraction 


interval 


of galaxies 


of galaxies 


0.00 - 0.45 




808 


14.54 % 


0.45 - 0.81 




998 


17.96 % 


0.81 - 1.11 




885 


15.92 % 


1.11 - 1.61 




898 


16.16 % 


1.61 - 2.15 




504 


9.07 % 


2.15 - 2.91 




746 


13.42 % 


2.91 - 4.01 




549 


9.88 % 


4.01 - 5.01 




150 


2.70 % 


> 5.01 




18 


0.32 % 


unknown 




2 


0.04 % 



also show the uncorrected luminosity function in the vari- 
ous plots as open circles) . The redshift binning was chosen 
such that we have good statistics in every redshift bin and 
that the influence of redshift clustering was minimized. 
The redshift binning and the number of galaxies in every 
bin is shown in Tabled 

The errors of the luminosity functions are calculated 
by means of Monte-Carlo simulations as follows. The pho- 



tometric redshift code provides redshift probability distri- 
butions P(z) for each single galaxy. In each Monte-Carlo 
realization, we randomly pick a new redshift for each ob- 
ject from a sample of redshifts distributed like P(z) and 
calculate the corresponding luminosity. This we repeat 250 
times which allows us to derive the dispersion of the galaxy 
number density z) for each magnitude and redshift 

bin due to the finite width of P(z) for each galaxy. The 
total error in (j) is finally obtained by adding in quadra- 
ture the error from the Monte-Carlo simulations and the 
Poissonian error derived from the number of objects in the 
bin. 

Photometric redshift errors may, in principle, affect the 
shape of the luminosity function at the bright end: By 
scattering objects to higher redshifts they let the steep 
fall-o ff at high luminosities appear shallower toro rv et all 
|2003() . However, in the case of the FDF the redshift errors 
are so small that the influence on the shape of the lumi- 
nosity function is negligible. 

5.2. The slope of the luminosity function 

We first investigate the redshift evolution of the faint-end 
slope of the luminosity function by fitting all three param- 
eters of the Schechter function (M*, (p* , and a). The best 
fitting a and the corresponding Icr errors for all wavebands 
and redshifts are listed in Table [3 

Despite of the depth of the FDF, Table[51shows that it 
is only possible to obtain reasonably tight constraints on 
the slope a for z < 1.5. In addition, strong parameter cou- 
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Fig. 4. Left panel: Absolute B magnitudes of galaxies in the FDF against redshift. Colors indicate spectral types 
(red to blue: old to young). The two lines indicate the 50% completeness limit for a red and a blue spectral energy 
distribution corresponding to an I-band limiting magnitude of 26.8. Right panel: Redshift number distribution of all 
galaxies in the FDF sample. The clustering observed in photometric redshift space (both panels) is probably mostly 
real, as we see clustered spectroscopic redshifts at z = 0.22, z = 0.33, z ~ 0.39, z = 0.45, z = 0.77, z = 2.35 and 
possibly at z = 0.95, z = 3.15, and z = 3.4. 



Table 2. Slope of the luminosity function for all wavelengths and all redshifts as derived from 3-parameter Schechter 
fits. 



z 


a (1500 A) 


Q (2800 A) 


a (u') 


" (g') 


a (B) 


[0.45,0.81] 


-1.14 (+0.08 


-0.07) 


-1.23 (+0.08 


-0.07) 


-1.27 (+0.06 


-0.05) 


-1.34 


(+0.05 


-0.03) 


-1.30 (+0.05 


-0.03) 


[0.81,1.11] 


-0.96 (+0.13 


-0.10) 


-0.99 (+0.10 


-0.08) 


-0.93 (+0.09 


-0.07) 


-1.16 


(+0.07 


-0.04) 


-1.21 (+0.07 


-0.04) 


[1.11,1.61] 


-1.05 (+0.18 


-0.16) 


-1.03 (+0.13 


-0.11) 


-0.95 (+0.10 


-0.09) 


-1.13 


(+0.11 


-0.09) 


-1.12 (+0.09 


-0.07) 


[1.61,2.15] 


-0.81 (+0.48 


-0.45) 


-0.97 (+0.32 


-0.28) 


-0.80 (+0.31 


-0.27) 


-1.29 


(+0.24 


-0.21) 


-1.33 (+0.27 


-0.20) 


[2.15,2.91] 


-0.38 (+0.21 


-0.15) 


-0.67 (+0.18 


-0.15) 


-0.70 (+0.16 


-0.16) 


-0.89 


(+0.22 


-0.15) 


-0.70 (+0.24 


-0.21) 


[2.91,4.01] 


-0.98 (+0.28 


-0.24) 


-0.95 (+0.19 


-0.17) 


-1.25 (+0.19 


-0.14) 


-1.24 


(+0.23 


-0.20) 


-1.30 (+0.27 


-0.20) 


[4.01,5.01] 


-0.77 (+0.38 


-0.26) 


-1.03 (+0.46 


-0.35) 


-1.09 (+0.54 


-0.27) 


-1.18 


(+0.37 


-0.21) 


-0.77 (+0.49 


-0.39) 



pling is observed between M* and a (see Fig. IC.ll in the 
Appendix[nj . We find only marginal evidence for a change 
of a with redshift for all wavebands. The lowest redshift 
bin (0.15 < z < 0.45), which we excluded from the fit 
because of poor number statistics in bright objects, gen- 
erally shows the steepest faint-end slope. Beyond redshift 
0.5, all data are consistent with a constant and shallow 
faint-end slope. 

We obtain as best error-weighted values for all red- 
shifts between 0.45 and 5.0 the numbers given in Table 01 
(upper part), assuming that a does not depend on red- 
shift. The slopes in the 1500 A, 2800 A, and u' band are 
very similar. The same applies for the slope in the g' and 
B band. Therefore, we combined the data for the 1500 A, 
2800 A, and u' band as well as for the g' and B band 
and derived combined slopes with an error weighted fit to 



the data of Table El The results are also listed in Table O 
(lower part). 

Almost all of the slopes listed in Table (21 are compat- 
ible within 2(7 with the slopes in Table (31 Therefore, we 
fixed the slope to these values for further analysis. This 
simplification is also justified by the fact that for all sub- 
sequent fits with fixed slope the reduced was generally 
close to 1. 

As a last test, we investigated the influence of the red- 
shift binning on the slope a. We enlarged our first two 
redshift bins to 0.41 < z < 0.91 (1433 galaxies) and 
0.91 < z < 1.61 (1438 galaxies) which allowed us to deter- 
mine luminosity functions with lower errors in all wave- 
bands. The slopes derived in these two larger bins were 
compatible with our previously derived fixed slope in ev- 
ery waveband. 
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Fig. 5. Luminosity functions at 2800 A from low redshift ((z) — 0.3, upper left panel) to high rcdshift {{z) = 5.5, 
lower right panel). The filled (open) symbols show the luminosity function corrected (uncorrected) for V/Vmax- The 
fitted Schechter functions for a fixed slope a are shown as solid lines. Note that we only fit the luminosity functions 
from (z) = 0.6 to (z) — 4.5. The parameters of the Schechter functions are given in 

Table E21 The Schechter fit for 

redshift {z) = 0.6 is indicated as dashed line in all panels. 



5.3. The restframe luminosity functions 

In this section we analyze the luminosity function in the 
UV (1500 A, 2800 A), u', g', and B band by means of a 
Schechter function fit with fixed slope (see Sect. I5.2|l . 

In the UV, we evaluate the luminosity function in 
two rectangular filters centered at 1500 ± 100 A and 
2800 ± 100 A. There are three reasons to analyze both 



wavelengths. First, for our lowest redshift bin ((z) ~ 0.6) 
the restframe magnitude derived at 2800 A is more robust 
than the one at 1500 A because the restframe wavelength 
of 2800 A corresponds to the observed U and does not need 
extrapolation to shorter wavelength. Second, we also in- 
clude the 1500 A luminosity function as it corresponds to 
a frequently used reference wavelength and is very well de- 
termined beyond redshifts of 2.5. Third, we want to show 
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Fig. 6. Luminosity functions in the g'-band from low redshift ((z) = 0.3, upper left panel) to high redshift ((z) = 5.5, 
lower right panel). The filled (open) symbols show the luminosity function corrected (uncorrected) for V/Vmax- The 
fitted Schechter functions for a fixed slope a are shown as solid lines. Note that we only fit the luminosity functions 
from (z) — 0.6 to (z) — 4.5. The parameters of the Schechter functions can be found in Tabl e The dotted line 
represents the local g'-band luminosity function derived from the SDSS ijBlanton et The Schechter fit for 

redshift (z) = 0.6 is indicated as dashed line in all panels. 



that the galaxy luminosity functions at both wavelengths 
are very similar and show the same redshift evolution. 

In the optical bands, we calculated the evolution of the 
lu minosity functions in the u' and g' bands (g' of SDSS, 
see lFukugita et al1ll996l not to be confused with Gunn g 
which was part of the filter set with which we observed 
the FDF). Because many authors have already published 



luminosity functions in the Johnson B-band, we include 
also this filter in our analysis. 

In Figs.|5land|Hlwe present the luminosity functions at 
2800 A and in the g' band, while the results at 1500 A as 
well as for the u' and B bands can be found in Figs. lA.ll 
IA.2l and lA.3l in appendix 1X1 The filled (open) symbols de- 
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Fig. 7. Icr and 2a confidence levels in Schechter parameter space for the different redshift bins, labeled by their mean 
redshift. A Schechter function with a fixed slope (see Table |2J has been fitted to the luminosity function at 1500 A 
{upper left panel), 2800 A {upper right panel), in the u'-band {lower left panel) and in the g'-band {lower right panel). 
The parameters of the Schechter function can be found in Table EH Table fOl Table FOl and Table EH 



note the luminosity function with (without) completeness 
correction. 

Even without fitting Schechter fimctions to the data, 
it is obvious that there is strong evolution in characteristic 
luminosity and number density between redshifts 0.6 and 
4.5. 

The solid lines show the Schechter function fitted to 
the luminosity function. The best fitting Schechter param- 
eter, the redshift binning as well as the reduced are also 
listed. The reduced ai'e quite good for all but one red- 
shift bin (2.15 < z < 2.91). The slope we adopted is not 
suitable for that bin and increases the x^. The depth of 
the FDF allows us to trace the luminosity function over a 
relatively large magnitude range. Even in our highest red- 
shift bin (4.01 < z < 5.01) the luminosity function spans 
an interval of 4 magnitudes. 



In Fig.[7|we show the la and 2a confidence contours of 
M* and 4>* for the different redshift bins, illustrating the 
correlation of the two Schechter parameters. The contours 
correspond to Ax^ — 2.30 and Ax^ = 6.17 above the 
minimum ■ The best fitting Schechter parameters and 
their la errors are summarized in Tables lA.ll IA.2I IA.3I 
El and 1131 for the 1500 A, 2800 A, u', g' and B bands, 
respectively. The la errorbars of the single parameters 
are derived from the projections of the two-dimensional 
contours using Ax^ — 1- 

We find a systematic brightening of M* and a system- 
atic decrease of cj)* from low to high redshift. The evolu- 
tion is very strong at 1500 A (upper left panel), 2800 A 
(upper right panel) and in the u'-band (lower left panel) 
and moderately strong in the g'-band (lower right panel). 
We do not show the B-band results as they behave almost 
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Table 3. Slope a of the luminosity functions for the dif- 
ferent wavebands as determined from an error-weighted 
fit to the data in Table [21 under the assumption that 
a{z) = const, (upper part). In the lower part of the Table 
we show the best values of a after combining the UV bands 
and the blue optical bands. 



filter 


a(z) = const 


1500 A 


-1.01 ±0.08 


2800 A 


-1.06 ±0.07 


u' 


-1.10 ±0.08 


g' 


-1.26 ±0.04 


B 


-1.24 ±0.04 


1500 A & 2800 A & u' 
g' &B 


-1.07 ±0.04 
-1.25 ±0.03 



identical as the g'-band. Although the variation of M* and 
(jf between adjacent redshift bins is in part influenced by 
large scale structure, the overall trend in the evolution of 
M* and cff is very robust. 

Since the integral of the luminosity function in the 
UV is strongly rela ted to the star-formation rate (SFR) 
l|Madaii et al.Lll998[l . we can derive the star- formation his- 
tory from the evolution of the luminosity function. The 
brightening of M* and decrease of 0* in the UV leads to 
an increase of the SFR between 0.5 < z < 1.5, whereas it 
remains approximately constant between 1.5 < z < 4.0. 
A detailed analysis of the star-formation history will be 
presented in a future paper (Gabasch et al., in prepara- 
tion) , preliminary results are published in iGabasch et alJ 

iool. 



6. Parameterizing the evolution of the luminosity 
function 

In order to quantify the redshift evolution of M* and 
we assume the simple relations of the form: 



M*iz) ^ M* + aln(l±z), 
0* (z) = 0* (1 + z)'' , and 
a (z) = ao = const. 



(1) 



Parameterizing M*(z) = Mq + aln(l + z) is equiv- 
alent to a dependence of L* (z) = Lq{1 + z)^ with 
^ = -0.41n(10)a « -0.921a. 

The best fitting values for a, 6, Mg , and c/jq are derived 
by minimizing 

= xHa,b,M*,^*) 

^ [m^J)~'^{M^,,z„a,b,MS,<j>*o)f ^ ^2) 

j=l i=\ 
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Fig. 8. 1(7 and 2a confidence levels of the parameters a 
and b in diflferent bands (1500 A, 2800 A, u', g' and B) for 
the evolutionary model described in the text. The values 
for a and 6 can be found in Tabled 



for the galaxy number densities in all magnitude 
and redshift bins simultaneously. 4>{Mij) is the number 
density of galaxies in magnitude bin i at redshift Zj\ 
'i'{Mij, Zj,a,b,MQ,(j)l) is the Schechter function evolved 
to the redshift Zj according to the evolution model defined 
in equation Q), and ct^ is the rms error of the luminosity 
function value. The resulting values for a, b, Mq, and (f>Q 
can be found in Table 0] 

The la and 2a confidence levels of the evolution pa- 
rameters a and b are shown for the different filters in 
Fig.lHl These contours were derived by projecting the four- 
dimensional distribution to the a-b plane, i.e. for given 
a and b we use the value of AIq and (/)q which minimizes 
the x^io-i b)- 

In Fig. 1^1 we show the relative redshift evolution of 
M* (left panel) and 0* (right panel) in the chosen filters. 
The Schechter parameters are the ones given in the tables 
in Appendix ^ The solid lines show the relative change 
according to our evolutionary model with the parameters 
from Table H 

Note that a, 6, Mq, and 0q were derived by minimizing 
equation (0) and not the differences between the (best 
fitting) lines and the data points in Fig. |51 

Fig.|51shows that the simple parameterization we have 
chosen with equation ^ describes the evolution of the 
galaxy luminosity functions very well. Still, the reduced 
xi values are somewhat larger than unity ('^ 4), because 
our approximations for evolution and faint-end slope may 
not be adequate for every redshift bin and because of the 
influence of large scale structure. Nevertheless, as there 



12 



Gabasch et al.: The evolution of the luminosity functions in the PDF: I. The blue bands 



Table 4. Evolution parameters according to equation 



filter 



1500 A 
2800 A 



B 



MS 
(mag) 



(10~^Mpc" 



-2.19 
-2.05 
-1.80 
-1.08 
-1.03 



-1-0.19 
0.19 

-F0.23 
0.24 

+0.24 
-0.21 

-1-0.30 
0.28 

+0.23 
0.28 



-1.76 
-1.74 
-1.70 
-1.29 
-1.27 



+0.15 

0.15 

+0.15 
0.16 

+0.14 
0.15 

+0.18 
0.18 

+0.16 
0.19 



-17 40+°-^^ 



9 71 +0.47 
^- ' ^-0.38 



+0.27 
0.26 

+0.24 
0.26 



-21.001°:^,? 



-18.16 
-18.95 



2.46: 
2.19" 



0.37 
0.37 



0.83trJ 



-20.92 



+0.32 
0.25 



0.82 



+0.14 
0.12 



are no stringent theoretical predictions for the evolution 
of M* and 0* we do not want to increase the number of 
free parameters, but increase the errors of a, 6, Mq, and 
(/)q instead. We do this by an appropriate scaling of the 
errors Uij of equation to obtain a of unity. 

For comparison, we also show in Fig. El the local val- 
ues from the SDSS ^Blanton et al.U200lj) . There is good 
agreement in the u'-band for both M* and 0* between our 
extrapolated values and the SDSS values. In the g'-band 
the value of M* is lower than the predicted one, but still 
within the Icr error of the Mn . 



7. Comparison with the literature 

In this section we compare the luminosity functions 
derived in the FDF with the luminosity functions of 
other surveys. As the cosmology and the wavebands in 
which the luminosity functions were determined are dif- 
ferent from ours for most of the surveys we chose the 
following approach. First we convert results from the 
literature to our cosmology {^m — 0.3, ilA = 0.7 and 
TJq = 70 kms~^ Mpc~^). Note that this conversion may 
not be perfect, because we can only transform num- 
ber densities and magnitudes but lack the knowledge of 
the individual magnitudes and redshifts of the galaxies. 
Nevertheless, the errors introduced in this way are not 
large and the method is suitable for our purpose. Second^ 
in order to avoid uncertainties due to conversion between 
different filter bands, we always use the same band as the 
survey we want to compare with. Thirds we also try to use 
the same redshift binning if possible. In addition, if the 
number of galaxies in the FDF is too small to derive a 
well sampled luminosity function we increase the binning. 

To visualize the errors of the literature luminosity func- 
tions we perform Monte-Carlo simulations using the AM*, 
A0*, and Aa given in the papers. In cases where not all of 
these values could be found in the paper, this is mentioned 
in the figure caption. We do not take into account any cor- 
relation between the Schechter parameters and assume a 
Gaussian distribution of the errors AM*, Ac/)*, and Aa. 
From 1000 simulated Schechter functions we derive the 
region where 68.8 % of the realizations lie. The resulting 
region, roughly corresponding to Icr errors, is shaded in 
the figures. The luminosity functions derived in the FDF 



are also shown as filled and open circles. The filled circles 
are completeness corrected whereas the open circles are 
not corrected. The redshift binning used to derive the lu- 
minosity function in the FDF is given in the lower right 
part of every figure. Moreover, the limiting magnitude of 
the respective survey is indicated by the low-luminosity 
cut-off of the shaded region in all figures. If the limiting 
magnitude was not explicitly given it was estimated from 
the figures in the literature. 

We first co mpare our lurninosit y functions in the UV 
to the results of Steidel et al.l lll999l) and the spectroscopic 
studies of Wi lson et all l)2002|) . 

Fig. ^1 shows a co mparison of th e 1700 A luminos- 
ity function derived by ISteidel et alJ (^99) at redshift 
(z) ~ 3.04 (left panel) and (z) ^ 4.13 (middle panel) with 
the luminosity func tion in the FDF. The galaxy sample of 
ISteidel et all l)l999|) is based on a R-band ((z) ~ 3.04) and 
an I-band ((z) ^ 4.13) selected catalogue and therefore 
similar to our I-band selected sample. Candidate galaxies 
were identified with the Lyman-break technique and most 
of them spectroscopically confirmed (564 galaxies of the 
(z) ~ 3.04 and 46 of the (z) ^ 4.13 sample, respectively). 

To derive the associated erro rs (shaded reg i on) of the 
Schechter functions derived by ISteidel et al.l (ll999^ we 
use the errors of M* and a of the (z) ~ 3.04 sample as 
given in Fig. 8 of their paper. As there are no errors re- 
ported for the (z) ~ 4.13 sample we assume the same er- 
rors as for the (z) ~ 3.04 sample. Therefore, the shaded 
region in Fig. 1101 (middle panel) is probably a lower limit. 
For the luminosity function in the FDF we use a red- 
shift binning of 2.54 < z < 3.54 (789 galaxies), and 
3.70 < z < 4.56 (144 galaxies) with the mean redshift 
of (z) ^ 3.04 and ( z) ~ 4.13 to be as close as possible to 
ISteidel et alJ l)l999l) 's mean redshifts. 

Fig. EH (left and middle panel) shows that there is 
very good agreement between the r esults derived in the 
FDF and the luminosity function of ISteidel et alJ l(l999^ 
if we focus only on the luminosity function brighter than 
the limiting magnitudes (shaded regions). On the other 
hand, because of the depth of the FDF we can trace the 
luminosity function 2 magnitudes deeper and therefore 
give better constraints on the slope of the Schechter 
function. We show in Fig. ^| (right panel) the Icr and 
2cr confidence levels for M* and a for a 3 parameter 
Schechter fit as derived from the FDF in the redshift 
interval 2.54 < z < 3.54 (solid fine) and 3.70 < z < 4.56 
(dotted line) . The steep slope a — —1.6 derived by 
ISteidel et alJ 11999") as marked by the horizontal dashed 
line can be excluded on a 2cr level. 



IWilson et al.l l|2002f) used galaxies selected in the 
restframe UV with spectroscopic redshifts to derive 
the luminosity function at 2500 A in 3 redshift bins: 
0.2 < z < 0.5 (U'-selected; 403 galaxies), 0.6 < z < 1.0 
(B-selected; 414 galaxies) and 1.0 < z < 1.5 (V-selected; 
518 galaxies). As the sample is not deep e nough to con- 
strain the slope of the Schechter function IWilson et 
((20021 used two fixed slopes of a = —1.0 and a = —1.5 
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Fig. 9. Redshift evolution of M* (left panel) and cj)* (right panel) for the filters g' (filled squares), u' (open triangles) 
and the t wo UV bands at 2800 A and 1500 A (filled circles). The arrows mark the values for M* and (p* as derived in 
the SDSS l)Blanton et all 120011) . 




Fig. 10. Compari son of the luminosity function at 1700 A of the FDF with the Schechter function derived in 
\Steidel et all figgd) : (z) ^ 3.04 (left panel) and (z) ^ 4.13 (middle panel). The filled (open) symbols show the lu- 
minosity function corrected (uncorrected) for V/Vmax- The shaded region in all plots is based only on AM*, and 
Aa (a detailed discussion concerning the errors AM*, and Aa can be found in Sect. Ej), where the cut-off at low 
luminosity indicates the limiting magnitude of the sample. In the right panel the la and 2(t confidence levels for M* 
and a for a 3 parameter Schechter fit as derived in the FDF in the redshift interval 2.54 < z < 3.54 (solid contours) 
and 3.70 < z < 4.56 (dotted contours) are shown. The horizontal dashed line marks the slope a = —1.6 as derived in 
ISteidelet all ill 9991) . 



to derive the best-fitting Schechter parameters. Since 
the errors of those parameters are not reported in the 
paper we can only make qualitative statements about the 
consistency of their and our luminosity functions: Fig. 1111 
shows that in the low and intermediate redshift bin there 
is reasonable agreement with our data, whi le in contrast 
to our result, the Schechter functions of I Wilson et alJ 



("2002^ do not show a significant brightening of M* 
their highest redshift bin. 



Comparison of the FDF luminosity function with 
the Schechter functions derived in ISullivan et all 



IW olf et al. 
IPn1' '^t ali 0101^, 



?2m ^. lKa,shikawaet a,].' 7200: 
Jwa,ta et abl (|2003l). lOuchi et 
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Fi g. 11. Comparison of the luminosity function at 2500 A of the FDF with the Schechter function derived 
in \ Wilson etaU hOOi) : 0.2 < z < 0.5 (left panel), 0.6 < z < 1.0 (middle panel), and 1.0 < z < 1.5 (right panel). 
IWilson et al.1 l|2n02^ fixed the slope to a = —1.0 (thick hne) and a — —1.5 (thin line) and used only M* and 0* 
as free parameters to determine the Schechter functions. 



ll2003bl). iBlanton et al.l l|200l|) . iBlanton et aD l|2003() . 
and IPoH et alJ l)2003j) are presented in Appendix 
In general, we find good agreement at the bright end, 
where literature datasets are complete. Differences in 
the faint-end slope in some cases can be attributed to 
the shallower limiting magnitudes of most of the other 
surveys. 

8. Comparison with model predictions 

As discussed in Sect. ^ key physical processes are 
involved in shaping the bright and the faint-end of the 
galaxy luminosity function. Therefore, it is interesting to 
compare luminosity functions predicted by models with 
observational results to better constrain those processes. 
In this section we compare the B-band luminosity func- 
tion in different r edshif t bi ns with mode l pred ictions of 
iKauffmann et all and lMenci et alJ l)2002() . 

\KaufFmann et al\ fiffPil) ; 

In Fig. El we show the B-band luminosity function of the 
F DF together with t he se mi-analytic model predictions 
bv IKauffmann et al.l l)l999(l ^ for the following redshifts: 
(z) ~ 0.20, (z) ~ 0.62, (z) - 1.05, (z) - 1.46, (z) - 2.12, 
and (z) ~ 2.97. 

There seems to be reasonably good agreement between 
the models (solid lines) and the luminosity functions 
derived in the FDF up to redshift (z) ^ 2.12. Of course at 
z « the model is tuned to reproduce the data. At z 3, 
the discrepancy increases as the model does not contain 
enough bright galaxies. Unfortunately, the models only 
predict luminosities for massive galaxies and, therefore, 
they do not predict galaxy number densities below M* . 



^ The models were taken from: 
http:/ /www.mpa-garching.mpg.de/Virgo/data_download.html| 



\Menci et al] \200i) : 

In Fig. 1131 we compare the B-band luminos ity functions 
of the FDF with the semi-analytic model bv lMenci et alJ 
l|2002|) for the following redshifts: (z) - 0.3, (z) ~ 0.6, 
(z) ~ 0.9, (z) - 1.4, (z) - 1.9, (z) - 2.6, (z) - 3.4, and 
(z) ~ 4.3. 

The agreement between the FDF data and the model 
in the lowest redshift bin (z) ~ 0.3 is very good, but this 
is again expected (see comment above). Moreover, if one 
focuses the co mparison only on t he hig her luminosity bins 
considered bv lKauffmann et al.l l)l999 V similar acceptable 
agreement with the data is observed. However, at lower 
luminosities and higher redshifts, the galaxy density of 
the simulation is much higher than the observed one. 



9. Summary and conclusions 

We analyzed a sample of about 5600 I-band selected galax- 
ies in the FORS Deep Field down to a limiting magnitude 
of I = 26.8 mag. 

A com parison with the very deep K-selected catalogue 
of lLabbe~et al.. L2003,) shows, that more than 90 % of their 
objects are brighter than our limiting I-band magnitude. 
Therefore our scientific conclusions arc not affected by this 
color bias. 

Based on 9 filters we derived accurate photometric 
redshifts with Az/(z.q ppp -I- 1) w 0.03 if compared with the 
spectroscopic sample ||No11 et aLi..2003.; Bohm et al.. 2003|) 
of 362 objects. We calculated and presented the luminos- 
ity functions in the UV (1500 A and 2800 A), u', B, and 
g' bands in the redshift range 0.5 < z < 5.0. The error 
budget of the luminosity functions includes both the pho- 
tometric redshift error as well as the Poissonian error. 

The faint-end slope of the luminosity function docs not 
show a large redshift evolution and is compatible within 
2a with a constant slope in most of the redshift bins and 
wavelengths considered here. Furthermore, the slope in 
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Fig. 1 2. Comparison of the B-band luminosity function of the FDF with predictions based on \Kauffmann et ai\ 
\l99^) {so\iA line): (z) - 0.20, {z) - 0.62, {z) - 1.05, (z) - 1.46, (z) - 2.12, and (z) ~ 2.97 (from upper left to lower 
right panel). The filled (open) symbols show the luminosity function corrected (uncorrected) for V/ Vmax.- The drops of 
the theoretical curves towards the faint-end is caused by the limited mass resolution of the models, see lKaufFmann et alJ 
1,1999.) for details. 



the 1500 A, 2800 A, and u' band is very similar but differs 
from the slope in the g' and B band. We derive a best 
fitting slope of a = — 1.07± 0.04 for the combined 1500 
A, 2800 A and u' bands and a = -1.25 ± 0.03 for the 
combined g' and B bands. We find no evidence for a very 
steep slope (a < — 1.6) at z ~ 3 and 1700 A rest wave- 
length as reported by o ther authors (e.g., ISteidel et al.l 
ll999ilOuchi et alJl2003bt) . From our data we can exclude 
a slope of a < —1.6 at redshift (z) ~ 3.0 and (z) ~ 4.0 at 
the 2(7 level. 

We investigate the evolution of M* and </>* by 
means of a redshift parameterization of the form 
M*(z) = Mq* + aln(l + z) and 0*(z) = 05(1 + z)^ We 
find a substantial brightening of M* and a decrease of 
with redshift in all analyzed wavelengths. If we follow the 
evolution of the characteristic luminosity from (z) ^ 0.5 
to (z) ~ 5, we find an increase of ^ 3.1 magnitudes in the 
UV, of ~ 2.6 magnitudes in the u' and of 1.6 magnitudes 
in the g' and B band. Simultaneously the characteristic 
density decreases by about 80 % - 90 % in all analyzed 
wavebands. 



Moreover, we compare the luminosity function derived 
in the FDF with previous observational datasets, mostly 
based on photometric results, and discuss discrepancies. In 
general, we find good agreement at the bright end, where 
their samples are complete. Differences in the faint-end 
slope in some cases can be attributed to the shallower lim- 
iting magnitudes of most of the other surveys. The only 
observations which reach the same limiting magnitudes as 
the FDF observations are th ose of Poh ct al. (2001 200311 
and the K-selected sample o f lKashikawa et al.l ' l)2003(l . The 
FDF results for the faint-end slope are in excellent agree- 
ment with those of Kashikawa et al. (2 00^1) but th e slope 
of th e Schechter function favored by IPoIT et al. 1 1I2OOII 
I2OO3I) is steeper than we would expect from the FDF. 

The semi-analytical models predict luminosity func- 
tions which describe (by construction) the data at low 
redshift quite well, but show growing disagreement with 
increasing redshifts. 
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results. Moreover, we thank Dr. N. Menci for providing an 
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Fig. 13. Compari son of the B-band luminosity function of the FDF with predictions based on the CDM model of 
\Menci et aD ^200^) (private communication; solid line): (z) 0.3, (z) - 0.6, (z) - 0.9, (z) ~ 1.4, (z) ~ 1.9, (z) ~ 2.6, 
(z) 3.4, and (z) ~ 4.3 (from upper left to lower right panel). The filled (open) symbols show the luminosity function 
corrected (uncorrected) for V/Vmax- 
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Fig. A.l. Luminosity functions at 1500 A from low redshift ((z) = 0.3, upper left panel) to high redshift ((z) — 5.5, 
lower right panel). The filled (open) symbols show the luminosity function corrected (uncorrected) for V/Vmax- The 
fitted Schechter functions for a fixed slope a are shown as solid lines. Note that we only fit the luminosity functions 
from (z) = 0.6 to {z) = 4.5. The parameters of the Schechter functions can be found in Table ED The Schechter fit 
for redshift (z) = 0.6 is indicated as dashed line in all panels. 
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Table A.l. Schechter function fit at 1500 A 



redshift interval 


M* (mag) 


ct>* (Mpc- 


') 


a (fixed) 


0.45 


- 0.81 


-18.17 +0.11 


-0.11 


0.0110 +0.0007 


-0.0006 


-1.07 


0.81 


- 1.11 


-18.85 +0.10 


-0.10 


0.0103 +0.0006 


-0.0006 


-1.07 


1.11 


- 1.61 


-19.48 +0.11 


-0.11 


0.0056 +0.0006 


-0.0005 


-1.07 


1.61 


- 2.15 


-19.97 +0.22 


-0.24 


0.0033 +0.0006 


-0.0006 


-1.07 


2.15 


- 2.91 


-20.61 +0.09 


-0.09 


0.0032 +0.0002 


-0.0002 


-1.07 


2.91 


- 4.01 


-20.72 +0.09 


-0.10 


0.0023 +0.0002 


-0.0002 


-1.07 


4.01 


- 5.01 


-21.00 +0.15 


-0.11 


0.0010 +0.0001 


-0.0001 


-1.07 



Table A.2. Schechter function fit at 2800 A 



redshift interval 


M* (mag) 


(j)* (Mpc- 


') 


a (fixed) 


0.45 


- 0.81 


-18.80 +0.15 


-0.15 


0.0104 +0.0007 


-0.0007 


-1.07 


0.81 


- 1.11 


-19.52 +0.09 


-0.10 


0.0089 +0.0005 


-0.0005 


-1.07 


1.11 


- 1.61 


-20.03 +0.09 


-0.09 


0.0053 +0.0004 


-0.0004 


-1.07 


1.61 


- 2.15 


-20.43 +0.18 


-0.17 


0.0029 +0.0005 


-0.0004 


-1.07 


2.15 


- 2.91 


-21.16 +0.09 


-0.08 


0.0030 +0.0002 


-0.0002 


-1.07 


2.91 


- 4.01 


-21.19 +0.10 


-0.08 


0.0021 +0.0002 


-0.0001 


-1.07 


4.01 


- 5.01 


-21.55 +0.17 


-0.21 


0.0009 +0.0001 


-0.0001 


-1.07 



Table A. 3. Schechter function fit in the u'-band 



redshift interval 


M* (mag) 


(Mpc- 


') 


a (fixed) 


0.45 


- 0.81 


-19.56 +0.16 


-0 


15 


0.0096 +0.0006 


-0.0006 


-1.07 


0.81 


- 1.11 


-20.12 +0.10 


-0 


10 


0.0080 +0.0004 


-0.0004 


-1.07 


1.11 


- 1.61 


-20.56 +0.08 


-0 


09 


0.0049 +0.0003 


-0.0003 


-1.07 


1.61 


- 2.15 


-20.70 +0.18 


-0 


16 


0.0033 +0.0004 


-0.0004 


-1.07 


2.15 


- 2.91 


-21.50 +0.08 


-0 


08 


0.0032 +0.0002 


-0.0002 


-1.07 


2.91 


- 4.01 


-21.57 +0.11 


-0 


10 


0.0022 +0.0002 


-0.0002 


-1.07 


4.01 


- 5.01 


-21.92 +0.24 


-0 


20 


0.0008 +0.0002 


-0.0001 


-1.07 



Table A. 4. Schechter function fit in the g'-band 



redshift 


interval 


M* (mag) 


(j)* (Mpc" 




a (fixed) 


0.45 - 


0.81 


-21.47 +0.20 


-0.20 


0.0042 +0.0003 


-0.0003 


-1.25 


0.81 - 


- 1.11 


-21.72 +0.15 


-0.15 


0.0039 +0.0003 


-0.0003 


-1.25 


1.11 - 


- 1.61 


-22.01 +0.14 


-0.14 


0.0026 +0.0002 


-0.0002 


-1.25 


1.61 - 


- 2.15 


-21.82 +0.20 


-0.20 


0.0020 +0.0004 


-0.0003 


-1.25 


2.15 - 


- 2.91 


-22.62 +0.13 


-0.10 


0.0020 +0.0002 


-0.0002 


-1.25 


2.91 - 


- 4.01 


-22.51 +0.13 


-0.14 


0.0016 +0.0002 


-0.0002 


-1.25 


4.01 - 


- 5.01 


-23.12 +0.22 


-0.23 


0.0005 +0.0001 


-0.0001 


-1.25 



Table A. 5. Schechter function fit in the B-band 



redshift interval 


M* (mag) 


<j>* (Mpc- 


') 


a (fixed) 


0.45 


- 0.81 


-21.28 +0.21 - 


-0.18 


0.0042 +0.0004 


-0.0003 


-1.25 


0.81 


- 1.11 


-21.57 +0.15 - 


-0.13 


0.0040 +0.0003 


-0.0002 


-1.25 


1.11 


- 1.61 


-21.91 +0.13 - 


-0.13 


0.0024 +0.0002 


-0.0002 


-1.25 


1.61 


- 2.15 


-21.65 +0.22 - 


-0.22 


0.0021 +0.0004 


-0.0004 


-1.25 


2.15 


- 2.91 


-22.44 +0.11 - 


-0.09 


0.0021 +0.0002 


-0.0002 


-1.25 


2.91 


- 4.01 


-22.44 +0.15 - 


-0.14 


0.0015 +0.0002 


-0.0002 


-1.25 


4.01 


- 5.01 


-22.81 +0.21 - 


-0.25 


0.0005 +0.0001 


-0.0001 


-1.25 
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Fig. A. 2. Luminosity functions in the u'-band from low redshift ((2) = 0.3, upper left panel) to high redshift ((z) — 5.5, 
lower right panel). The filled (open) symbols show the luminosity function corrected (uncorrected) for V/Vmax- The 
fitted Schechter functions for a fixed slope a are shown as solid lines. Note that we only fit the luminosity functions 
from (z) — 0.6 to (z) — 4.5. The parameters of the Schechter functions can be found in Tabl e r^~3l The dotted line 
represents the local u'-band luminosity function derived from the SDSS ijBlanton et 1il l200lh . The Schechter fit for 
redshift (z) = 0.6 is indicated as dashed line in all panels. 
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Fig. A. 3. Luminosity functions in the B-band from low redshift ((z) = 0.3, upper left panel) to high redshift ((z) = 5.5, 
lower right panel). The filled (open) symbols show the luminosity function corrected (uncorrected) for V/Vmax- The 
fitted Schechter functions for a fixed slope a are shown as solid lines. Note that we only fit the luminosity functions 
from (z) = 0.6 to (z) = 4.5. The parameters of the Schechter functions can be found in Table Ol The Schechter fit 
for redshift (z) = 0.6 is indicated as dashed line in all panels. 
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Appendix B: Comparison with literature 

In this appendix we compare the luminosity functions de- 
rived in the FDF with the results of further publications 
as introduced in Sect.[7| The filled (open) circles show the 
completeness corrected (uncorrected) luminosity function 
as derived in the FDF in the redshift bin listed in the lower 
right corner. The solid line(s) represent the Schechter func- 
tion given in the different papers transformed to our cos- 
mology. To visualize the errors associated to this Schechter 
function we perform a Monte-Carlo simulation using the 
errors of the Schechter parameters reported in the specific 
paper (see Sect. 13 for more details). As the errors for all 
three Schechter parameters (AM*, Ai^*, and Aa) are not 
always given in the paper, we denote in the caption the 
errors used to perform the simulation. The region wherein 
68.8 % of the realizations lie are shown as shaded region 
in the plots and corresponds roughly to the la error due 
to the Schechter errors reported in the figure captions. 
Moreover the cut-off of the shaded region marks the lim- 
iting magnitude of the survey we compare with. 

B.l. UV bands 

\Sullwan et nl\ f^M): 

Although the volume of the FDF at low redshift is rather 
small, and therefore is not well suited to properly sample 
the bright end of the Schechter function, we compare for 
completeness in Fig. IB. II our luminositv function also 
with the luminosity function derived in ISiillivan et alJ 
l)200nl) . Their sample contains 433 UV-selected sources 
within an area of 2.2deg^. 273 of these objects are 
galaxies and cover the redshift range z ~ — 0.4. The 
solid line in Fig . IB. 1| represents the luminosity function 
at 2000 A from ISullivan et all l)2000|) whereas the filled 
circles shows our V/Vmax corrected luminosity function 
derived at 0.15 < z < 0.4. Despite the small volume, the 
I-selected catalogue and the extrapolated 2000 A lumi- 
nosity function (see above) there is a general agreement 
with only small systematic offsets (probably also due to a 
known cluster at z - 0.33 |(No11 et all 120031) ). This is an 
additional confirmation of the validity of our technique to 
derive the luminosity function as described in Sect. 15. ll 



\Wolf etal\ fiM): 

In Fig. lB.2l we compare the luminosity function at 2800 A 
of the FDF with the R-band select ed luminosity func tion 
derived in the COMBO-17 survey llWolf et all l20n3^ for 
different redshift bins: 0.2 - 0.6, 0.6 - 0.8, 0.8 - 1.0, 1.0 - 
1.2. Because of the limited sample size of the FDF at low 
re dshift we could no t use the same local redshift binning 
as IWolf et all l)2003|) . We compare therefore in Fig. lR2l 
(upper left panel) the C0MB017 Schechter function at 
{z) ^ 0.3 (light gray) and (z) ^ 0.5 (dark gray) with the 
FDF luminosity function derived at 0.2 < z < 0.6. There 
is an overall good agreement between the FDF data and 




Fig. B.l. Comparison of the luminosity function at 
20 00 A of the FDF wi th the Schechter function derived 
m lSullivan eTd\ \200A ) (z ~ ~ 0.4). The shaded region 
is based on AM*, A(/)*, and Aa, where the cut-off at low 
luminosity indicates the limiting magnitude of the sample. 

the COMBO-17 survey at all redshifts under investiga- 
tions if we compare only the magnitude range in common 
to both surveys (shaded region). Nevertheless the number 
density of the FDF seems to be slightly higher which 
mos t probably can b e attributed to cosmic variance. 
The IWolf et ahl l)2003l) team derived the faint-end slope 
from relatively shallow data which have only a limited 
sensitivity for the faint-end slope. Thus, the di sagreement 
between the much deeper FDF data and the IWolf et all 
(^2003) results at z ~ 0.5 and for z > 1 does not come as 
a surprise. 

iKashikawa et 

In Fig. IB .31 we compare our luminosity function with 
the K-band sel ected 2000 A luminosity function of 
IKashikawa et al.l (^003) derived in the Subaru Deep 
Survey. They used photometric redshift to determine the 
distance for 439 field galaxies. There is a good overall 
agreement of the luminosity functions in the redshift 
bins 0.6 < z < 1.0, 1.0 < z < 1.5, 1.5 < z < 2.5. 
Only in the highest redshift bin (2.5 < z < 3.5) the 
number density derived in IKashikawa et ahP l|2003) is 
lower by a factor of about 2 when compared with the FDF. 

Poli et all \200i^ : 

Poli et alJ l|200l|) combined three pencil beam surveys as 
the HDFN, the HDFS and the Ne w Technology Telescope 
Deep Field l)Arnouts et al.l Il999j) reducing the influence 
of cosmic variance and derived the 1700 A luminosity 
function at 2.5 < z < 3.5. In Fig. IB.4I we compare the 
result with the luminosity function in the FDF. There is 
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Fig. B.2. Comp arison of the luminosity function at 2800 A of the FDF with the Schechter function derived in 
I Wolf et all \20nA) : 0.2 < z < 0.4 (upper left panel, light gray), 0.4 < z < 0.6 (upper left panel, dark grey), 0.6 < z < 0.8 
(upper right panel), 0.8 < z < 1.0 (lower left panel), 1.0 < z < 1.2 (lower right panel). The shaded regions of all plots 
are based on AM*, A(/)*, and Aa, where the cut-off at low luminosity indicates the limiting magnitude of the sample. 



very good agreement although the slope of the Schechter 
function [a = —1.37) is slightly steeper than we would 
expect from the FDF. 



Iwata et al 



Iwata et al. 



'200W: 

l|2003l ) analyzed about 300 galaxies in a 
575 square-arcmin field detected in the I and z band at 
redshift z 5, selected by means of the Lyman-break 
technique. They derived the luminosity function at 
1700 A statistically. We analize Table 3 of llwata et al.l 
l)2003|) with the same method as described in Sect. IS.ll to 
get approximate errors for M* and (f>* for a fixed slope 
of a = —1.5 (as given in the paper). From these AM* 



and A(j)* we calculate the shaded region of Fig. IB. 51 (left 
panel). Fig. IB. 51 (left panel) compares the luminosity 
function of Iwata et alJ l|2003j) with the luminosity func- 
tion of the FDF deri ved at 4.01 < z < 5.01. Although 
the number density of llwata et al.l l|2003l) at z 5 seems 
to be slightly lower than the number density derived in 
the FDF at (z) ~ 4.5 the overall agreement is rather 
good. On the other hand, part of this decrease in density 
may also be due to intrinsic evolution between redshift 
(z) ^ 4.5 and (z) ~ 5.0. According to our evolution model 
as derived in Sect. El we would expect a decrease of 0* of 
about 15 %. 
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Fig. B.3. Comparison of the luminosity function at 2000 A of the FDF with the Schechtcr function derived in 
\Kashikawa et nl\ \200^) : 0.6 < z < 1.0 (upper left panel), 1.0 < z < 1.5 (upper right panel), 1.5 < z < 2.5 (lower left 
panel), 2.5 < z < 3.5 (lower right panel). The shaded regions of all plots are based on AM*, A(/)*, and Aa, where the 
cut-off at low luminosity indicates the limiting magnitude of the sample. 



Ouchi et al 



Ouchietal 



MM): 



1 2003b() investigated photometric properties 



of about 2600 Lyman-break galaxies at z = 3.5 — 5.2. 
Based on this sample they derived the luminosity func- 
tion at 1700 A for three redshift bins: z = 4.0 ± 0.5, 
z = 4.7 ± 0.5, z = 4.9 ± 0.3. In Fig. EH (right panel) we 
compare their Schechter function for a fixed slope of a = 
—1.6 with the luminosity function of the FDF derived at 
4.01 < z < 5.01. The Schechter function for z = 4.0±0.5 is 
shaded in dark gray, the z = 4.7±0.5 Schechter Function is 
shaded light gray and the z — 4.9 ±0.3 Schechter Function 
is represented by the dashed line (n o errors repo r ted). It 
is difficult to compare the results of lOuchi et al. 1 ^l2003b^ 



with the FDF. Our data favor a less steep slope of the lu - 
minosity function than advocated bv lOuchi et ai]ll2003bl) . 



B.2. SDSS bands (u', g', ^-^u, ^-^g) 

In this section we want to compare the luminosity function 
in the FDF with the one from the SDSS. 

In Fig. ESI (left panel) and Fig. |R7| (left panel) we 
show the luminosity function derived in iBlanton et alJ 
(; 2001.') for z ~ 0.1 in the u' and g' band, respectively, as 
light shaded region. To make a more appropriate compari- 
son between our 'local' results derived at 0.15 < z < 0.45 , 
we evolve the Schechter function of .Blanton et al JliooJ) 
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Fig. B.4. Comp arison of the luminosity function at 1700 A of the FDF with the Schechter function derived in 
\Poli et al\ \200i) (2 .50 < z < 3.50 ). The shaded region is based only on AM*, and Aa, where the cut-off at low 
luminosity indicates the limiting magnitude of the sample. 
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Fi g. B.5. Left panel : Comparison of the luminosity function at 1700 A of the FDF with the Schechter function derived 
in Iwata et "dl 1200 A) (z ~ 5). The shaded region is based only on AM*, and A(/ )*. Right pan el: Comparison of the 
luminosity function at 1700 A of the FDF with the Schechter functions derived \n \Ouchi et al\ \2003l )'■ z = 4.0 ± 0.5 
(dark shaded), z = 4.7 ±0.5 (hght shaded), and z — 4.9 ±0.3 (not shaded; dashed line). Both shaded regions are based 
only on AM*, and A(/)*, where the cut-off at low luminosity indicates the limiting magnitude of the sample. 
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to (z) ~ 0.3 according to our evolutionary model de- 
scribed in Sect. We use for the u'-band the parameter 
a = —1.80 and h — —1.70 whereas for the g'-band we 
use a = —1.08 and h = —1.29. The evolved Schechter 
function is shown as dark shaded region in in Fig. IB. 61 
(left panel) and Fig. IB. 71 (left panel) for the u' and g' 
band, respectively. Despite the small volume of the FDF 
in the local redshift bin, the agreement is very good in 
both bands and especially in the g'-band. We therefore 
conclude that there is no hint of a possible systematic 
offset between the two datasets. 

In Fig. IB. 61 fright panel) and Fig. IB. 71 (right panel) we 
also sh ow the luminosity function derived in Bla nton et al.l 
l)2003(l for the blue-shifted filter o ^u and °-^g. Again, the 
light shaded region represents the (z) '^0.1 luminosity 
function whereas the dark shaded region shows the lu- 
minosity function evolved to (z) ~ 0.3. We use the same 
evolutio n parameter a s derived for u' and g'. The approach 
used by iBlanton et al . (2003) differs from those used in 
all other studies, inc luding ours and the previous SDSS 
l)Blanton et all l200lh results. It is therefore beyond the 
scope of the paper to explain the discrepancies. 



B.3. B-band 




Poll et al.l l|2003'l analyzed 1541 I-selected and 138 K- 
selected galaxies to construct the B-band luminosity func- 
tion up to redshift (z) ^ 3. A comparis on between the 
luminosity function of IPoli et a,l.l l)2003|) and the FDF 
is shown in Fig. IB. 81 for the redshift bins 0.4 < z < 0.7 
(upper left panel), 0.7 < z < 1.0 (upper right panel), 
1.3 < z < 2.5 (lower left panel) and 2.5 < z < 3.5 (lower 
right panel). 

In neither of the redshift bins an error for </)* is reported 
in the paper and therfore could not be included in the 
simulation of the shaded region. For the two lower redshift 
bins (0.4 < z < 0.7 and 0.7 < z < 1.0) the shaded region 
is based on AM* and A(/)* whereas in the high redshift 
bins (1.3 < z < 2.5 and 2.5 < z < 3.5) the shown error of 
the Schechter function (shaded region) is bas ed only on 
AM*. If this is taken into account, the results of lPoli et alJ 
l|2003l) are in good agreement with the FDF, but again, the 
slope of the Schechter function is too steep when compared 
with the FDF luminosity function. On the other hand the 
brightening of M* with redshift is present in both samples. 



Appendix C: Confidence levels for the slope 
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Fig. B.6. Left panel: Comparison of the u'-band luminosity function of the FDF with the Schechter function derive d in 
\BlqntQn et al. (2001 J at z ^ 0.1 (hght shaded). The dark shaded region shows the Schechter function of Blantpn et alJ 
Ij200ll) evolved according to our evolutionary model described in Sect. to redshift z ~ 0.3. The shaded regions 
are based on AM*, A0*, and Aa. Right panel: Comparison of the "'^u-band luminosity function of the FDF with 
the Schechter function d erived in ^Blant on et al. i 200o.) at z ~ 0.1 (light shaded). The dark shaded region shows 
the Schechter function of iBlanton et al.l (2003(l evolved according to our evolutionary model described in Sect. to 
redshift z ~ 0.3. The shaded regions are based on AM*, At/)*, and Aa, where the cut-off at low luminosity indicates 
the limiting magnitude of the sample. 




Fig. B.7. Left panel : Comparison of the g'-band luminosity function of the FDF with the Schechter functi on derived in 
\Blanton et all ^200 A) at z ~ 0.1 (light shaded). The dark shaded region shows the Schechter function of .Blanton et alJ 
evolved according to our evolutionary model described in Sect. to redshift z ~ 0.3. The shaded regions 
are based on AM*, A0*, and Aa . Right panel: Compa rison of the °-^g-band luminosity function of the FDF with 
the Schechter function d erived in \Blanton et al\ ^200^ at z ~ 0.1 (light shaded). The dark shaded region shows 
the Schechter function of iBlanton et al.l l|2003() evolved according to our evolutionary model described in Sect. to 
redshift z ~ 0.3. The shaded regions are based on AM*, A(f>* , and Aa, where the cut-off at low luminosity indicates 
the limiting magnitude of the sample. 
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Fig. B .8. Comparison of the B-band luminosity fmiction of the FDF with the Schechter function derived m \Poli et ol\ 
0.4 < z < 0.7 (upper left panel), 0.7 < z < 1.0 (upper right panel), 1.3 < z < 2.5 (lower left panel), and 
2.5 < z < 3.5 (lower right panel). The shaded region is based only on AM*, and Aa for 0.4 < z < 0.7, and 0.7 < z < 1.0, 
whereas for the 1.3 < z < 2.5, and 1.3 < z < 2.5 the shaded region is based only on AM*, where the cut-off at low 
luminosity indicates the limiting magnitude of the sample. 
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Fig. C.l. 1(7 and 2a confidence levels in Schechter parameter space. A Schechter function with three free parameters 
M*, 0*, and a has been fitted to the luminosity function at 1500 A {upper left panel), 2800 A (upper right panel), u' 
(lower left panel) and g'-band (lower right panel) and projected to the M* - a plane. The various contours in each 
panel correspond to the different redshift bins, ranging from (z) = 0.6 (low luminosity) to (z) — 3.5 (high luminosity). 
We alternate continuous and dotted lines for clarity. The dashed line marks the fixed slope (a(z) — const.) used to 
derive the luminosity functions in the different wavebands (see Table 01 lower part). 



